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^ ' Abstract 

^ . This paper is devoted to the study of eigen-sequences for some important operators 

acting on sequences. Using functional equations involving generating functions, we 
completely solve the problem of characterizing the fixed sequences for the Generalized 
Binomial operator. We give some applications to integer sequences. In particular we 
show how we can generate fixed sequences for Generalized Binomial and their relation 
' with the Worpitzky transform. We illustrate this fact with some interesting examples 

and identities, related to Fibonacci, Catalan, Motzkin and Euler numbers. Finally 
we find the eigen-sequences for the mutual compositions of the operators Interpolated 
' Invert, Generalized Binomial and Revert. 

1 Introduction 



> 
in 



In a previous paper, Barbero, Cerruti, and Murru [3] have studied the action of a gener- 
alization of the Binomial interpolated operator on linear recurrent sequences. Particular 
attention was given to the study of the fixed sequences for this operator, finding some inter- 
esting results about known sequences, like the Catalan numbers. As final remarks, we posed 
the question of studying the general problem of characterizing the sequences fixed by this 
operator. The research interest involving operators which act on sequences has an important 
aspect related to those sequences which are fixed, also called eigen-sequences. Bernstein and 
Sloane [4] have deeply studied fixed sequences, for the composition of several operators as 
Binomial, Invert and others. The Binomial inverse transformation has been especially stud- 
ied in this sense. For example. Sun [10] and Wang [11] have found interesting results about 
the fixed sequences for the Binomial inverse operator. Furthermore, Mattarei and Tauraso 
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[7] have studied several congruences using these invariant sequences. Since the Binomial and 
the Binomial inverse transform arise as a particular case of Generalized Binomial operator, 
the relevance of finding the corresponding eigen-sequences is clear. Here we give a complete 
answer for this question characterizing the generating functions of such fixed sequences. 

Definition 1. We define, over an integral domain R, the set 

S{R) = {a= {an)t=o ■ Vn a„ G R, ao = l] . 

Definition 2. For any element x E R, the Interpolated Invert operator I^^^ acts on a 
sequence a G S{R) giving the sequence J(^^(a) = b = {bn{x))^^Q, where the sequence b has 
generating function 

+ 00 

B{t) = Y,bnixr= -=\^ . (1) 

"=o 1-xtY, ant"- 

n=0 

Definition 3. For any element y G R, the Interpolated Binomial operator L*^^^ acts on a 
sequence a G S{R) giving the sequence L'^^\a) = I = ilniy))n=o^ where the n-th element of 
the sequence / is 

^n = i2(%'"'^r (2) 
The exponential generating function for the sequence / is 

n=0 ' n=0 j=0 ^ J J- J- 



being 



n=0 n=0 



so that (recalling that Oq = 1) we have the ordinary generating function 

L(t) = -A ^ ^ 



t \\-yt 
with A{t) = J2 «nt"+^ . 

n=0 

Definition 4. The Revert operator r] acts as follows: if we consider the sequences a = (a„)^^ 
G S{R) , b = {bn)n=o ^ ^i^) have 77(a) = & if and only if the sequences a and b satisfy 
the relations 

u = u{t) = Y: ant^+^ 

(3) 

t = t{u) = J2 bnU^^^ inverse series of u . 

n=0 
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Prodinger [8] defined a generalization of tlie Interpolated Binomial operator as 



L^'^^y\a) = b, b^ = J2('')hY-'ai (4) 

for any sequence a G S{R) and nonzero elements h,y ^ R. The definition can be extended 
when /i = or y = as follows: 

L(0'^)(a) = (?/"ao)+ro, L^'^'\a) = {h-aX=o, L^^'^H^) = (^o, 0,0,0, ...)• 
We observe that the sequence L^^'y\a) has exponential generating function 

Ah{t) = ey'A{ht) 

and ordinary generating function 

L.W = ( , (5) 



ht \l~yt 

as one can find with a little bit of calculation, where the functions A{t) and A{t) are the 
same introduced in the previous Definition 3. The authors [3] have characterized linear 
recurrent sequences of degree 2 fixed by this operator and this result could be extended to 
linear recurrent sequences of higher degree. However, as we pointed out in that article, there 
are fixed points for the operator L'-'^'^-' which are not linear recurrent sequences. In this paper 
we give a full explanation of this fact, via functional equations (see for example Aczel [2] 
and Kannappan [5]), and we find the fixed sequences related to the operators I^^^ o L^'^'y\ 
7] o L^y^^) and rj o /(^). 



2 Fixed sequences for Generalized Binomial operator 

We observe that a sequence a E S{R) , fixed by the operator L^^'y\ has an exponential 
generating function A{t) which must satisfy the functional equation 

A{t) = Ah{t) = e^'Aiht) (6) 

with the condition ^(0) = = 1 . The analysis of this equation allows us to state some 
results and conditions on parameters h and y, summarized in the following theorem. 

Theorem 5. When h ^ 1, h ^ —1, and y ^ the only sequence fixed by the operator L^^'y^ 

When h = 1 we trivially have fixed sequences if and only ify = 0, in other words the operator 
]^{h,y) i)ecomes the identical operator. When y = 0, h ^ 1, and h ^ —1 we have the only fixed 
sequence A = {oq = l,a„ = V?i>l}. When y = and h = —1 we have fixed sequences 
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A for any even exponential generating function A{t) with A{0) = 1. Finally, when y ^ 
and h = —1 the exponential generating function of a fixed sequence fulfills the relation 

where the function (f){u) satisfies the functional equation 

(/){u) = ucj) ^-^ with u> and 0(1) = 1 . (8) 

Proof. First of all we consider the case h ^ 1 , h ^ —1, and y 0. If we substitute e^* = u, 
or equivalently t = , we have 



Now we naturally define (j) = Ao g o log, where g(t) = ^, and equation (6) becomes 

(p{u) = u(f){u^) , (10) 

where the condition ^(0) = 1 corresponds to the condition 0(1) = 1 . Iterating the previous 
relation (10), we easily obtain 

0(u) = u(j){u'') = n^+V(u^') = u^-^^+^'(j){u''") = ■■■ = u^4){v!'") . 

Hence, if < 1, taking the limit when n — )■ +oo, we find 

1 1 

(j){u) = U^-h(j)[l) = u^-i^ , 

which implies that 

+00 / \ ri 

The case \h\ > 1 gives in a similar way the same result. When h = —1 and y ^ the 
functional equation (8) is straightforward from equation (10). If /i = 1 the equation (6) 
becomes 

Ait) = e^'Ait) , 

which is obviously satisfied by sequences in S{R) if and only if y = 0. Moreover, when y = 
and h = —1, clearly A{t) = A{—t), so the exponential generating function A{t) is an even 
function. Finally, when y = 0, h ^ 1, and h ^ —1 the functional equation 

A{t) = A{ht) , 

with A{0) = 1, has the only constant solution A{t) = 1, corresponding to the fixed sequence 
A = {ao = l,an = Vn > 1}. □ 

With the next theorem we characterize all the solutions G C°° (K'*') of (8) with 0(1) = 1 . 
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Theorem 6. The equation (8) has a solution cf) G (M"*") with 0(1) = 1 if and only if exists 
an R-valued function F eC°^ (M+ x M+) such that 

z;F(1,1) = 1 

a) F{u, v) = F{v, u) 

Hi) F{cu, cv) = cF{u, v) Vc G M"'". 

Then = F{u, 1) = F{l,u). 

Proof. Let us consider an M- valued function F G (M^ x M+). From the symmetry 
F{u, v) = F{v, u) and using the homogeneity of degree 1 of F{u, f ), we have 



uF[l,-\ =vF (l,- 
u/ \ V 



and consequently 



Fil^] = ^F(l^ 

u 



In particular, when we pose v = 1 the previous equality becomes 

1 



F{l,u) =uF\l 



u 





= -<P 






V 





and clearly (f){u) = F{l,u) = F{u,l) is a solution of equation (8) belonging to C°° (M^). 
Now if G C°° (M"*") is a solution of equation (8) we have 

V(ti, i;) G R+ X M+ (11) 

and we can define an M-valued function F G C°° (M^ x ]R+) by posing F(u, v) = vcj) (^) . 
This function satisfies i^(l, 1) = 1, is symmetric, in fact from relation (11) 

F{u,v) = v(f) ^— j = u(/) ^— j = F{v,u) , 

and homogeneous of degree 1 because 

F{cu, cv) = cv(j) ( — ) = cv(p ( — I = cF{u, v) Vc G M"*" . 
\cv / \v / 

□ 



3 Some remarks, examples, and applications to integer 
sequences 

All the previous results, with the same conditions on parameters y and h, can be easily 
generalized in order to find fixed sequences whose first element oq belongs to M*. We also 
observe that there are some straightforward methods to generate solutions of equation (6) 
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starting from known solutions, when h = —1 and the parameter y has an assigned value. A 
simple calculation proves that if the function A{t) satisfies equation (6) when h = —1 , also 
the product f{t)A{t), with an even function f{t), is a solution, and any linear combination 
of solutions is also a solution for the same assigned value of the parameter y. Finally if we 
have a set of functions Ak{t) satisfying Ak{t) = e^'=*^fc(— t), k = 1,2, ■■ ■ , n, then A{t) = 

n 

Y\ {Ak{t))°'^ , ak G N, satisfies A{t) = e '^^i A{—t) . At last we have to consider the 
fc=i 

problem related to the existence of fixed sequences starting with exactly k zeros, or in other 
words sequences with ^^(0) = when j = 0, 1, ■ ■ ■ , - 1 and ^('=)(0) 0. Differentiating 
/c-times both members of equation (6) with respect to the variable t, we have 



j=0 

and if A{t) is an exponential generating function of a fixed sequence starting with k zeros, 
evaluating this equality for t = 0, we obtain ^'^'^-'(O) = h''A^''\0). This relation is true when 
k is odd, if and only if /i = 1 and consequently y = 0, the trivial case, and when k is even 
we retrieve the non trivial case h = —1. So we have fixed sequences starting with k zeros 
for L^~^'y^ if and only if A; = 2m. This fact allows us to state that a fixed sequences with 2m 
starting zeros is a shift of 2m position of a fixed sequence with ao 7^ 0, whose exponential 
generating function is a solution of equation (6) when h = —1. This is a particular case 
of what we have pointed out before, because the exponential generating function of such a 
sequence is the product of a known solution of equation (6), when h = —1, with the even 
function f{t) = t^"^. 

After these remarks let us illustrate some straightforward application of Theorem 6 to integer 
sequences. If we consider any real number a and we take 



the sequence 



n=0 



2 

is fixed for L^~^'^\ In particular, when a = ^-y-^, the sequence becomes 

where L„ is the n-th Lucas number. Another interesting situation is a = |, in this case 

(f){u) = y/u 

corresponding to 

Ait) = . 
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Thus all the sequences having exponential generating functions 

Ait) = eP'l^{qt)t^'' g e M*, z/ G Z , (12) 
where the function I,,(x) is the modified Bessel function of the first kind (see Abramowitz 



,x)x 



and Stegun [1]) , are eigen-sequences for the operator L*^~^'^^\ because the functions /, 
are even functions. In the OEIS [9] there are a lot of sequences whose exponential generating 
functions correspond to one of those in equation (12) for suitable values oi p , q, and p. The 
Motzkin numbers A001006 and the central binomial coefficients A000984 have exponential 
generating functions 

A{t) = e*/i(2t)r^ and A{t) = e'^^Io{2t) 

and they are eigen-sequences respectively for the operators L^'^'"^^ and L^~^'^\ Moreover the 
sequences of the form 

(for example the central Delannoy numbers A001850 and A115865), where the polynomial 
Pn{x) is the n-th Legendre polynomial (see Abramowitz and Stegun [1]) , have exponential 
generating functions 

A{t) = e'=™*/o (k^/m{m - l)t) 
and they are invariant for the operators /,(-i.2fcm)_ 

As a final application we want to show how to generate many fixed sequences starting from 
the operator L^'^'y\ finding some interesting formulas. Let us consider an even sequence 

fc= (60,0,62,0,64,0,66,0,...) 

we know that the sequence a = L^^'^^b) has exponential generating function 

+~ 6„/i"r 



E 



n=0 



i.e., has an exponential generating function given by e^* multiplied by an even function. 
Thus, as immediate consequence of the Theorem 5, the sequence a is a fixed sequence for 
the operator L^~^''^y\ From Definition 4 we find a relation between the sequences a and 6 

n L«/2J / \ / • \ 

E E J 2 J (-l)'2-y--'^/.^'=62. = a., Vn > 1. (13) 
i=0 k=o ^ / ^ / 

This equality gives a connection with the Worpitzky transformation [6], which converts any 
sequence s into a polynomial sequence s{x) as follows: 

" fc / 1- \ 

W{s) = s{x), sn{x) = J2 E(-l)"' ( m ) '^'^^"^ + + 

fc=0 m=0 ^ ^ 
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We observe that the polynomial sequence s{x) is an Appell sequence, i.e., 

L^^'y\s{x)) = s{x + y). 

When we consider an even sequence b there is an unique sequence a, which we call the 
worpification of b such that 

W{a) = a{x), a{0) = b. 

Since the polynomial sequences obtained by the Worpitzky transformation are Appel se- 
quences, we have 

L^^^y\b) = L^^^y\a{0)) = aiy) 
and so we know that a{y) is fixed by L^~^''^y\ Let us consider the Catalan numbers, A000108 , 

C = (1, 1, 2, 5, 14, 42, 132, 429, 1430, . . .) . 

The sequence b of Catalan numbers interpolated with O's, i.e., 

b = (1,0,1,0,2,0,5,0,14,0,42,0,...) 

and the sequence M of Motzkin numbers A001006 satisfy the relation L^^'^\b) = M (see 
OEIS [9], A000108 , Paul Barry 2003). Thus Motzkin numbers are fixed by L^-^'^), as we 
have previously showed, and using (13) we have the following relation between Motzkin and 
Catalan numbers 

n L'/2J 



n 1,—n \ / \ / 



1=0 k=0 

In particular, using our techniques we have found that 

LV2J 



EGlf-*^- 

i—n N / 



k=0 

which is a known result and moreover we have 



^ ("'^^-1)^2"-*^ = M„. 



The last relation can be found directly as a consequence of the Corollary 3.3 of the work 
of Sun [10]. Furthermore, we have a connection between Catalan and Motzkin numbers 
in terms of the Worpitzky transformation. Indeed, if we call a the worpification of b and 
W{a) = a{x), then we have that 

a(0) = b, a(l) = M, a(2) = C. 

Another application involves, for example, the sequence F of Fibonacci numbers A000045 . 
The Binomial transform of the aerated Fibonacci numbers 

6 = (0, 0, 1, 0, 1, 0, 2, 0, 3, 0, 5, 0, 8, 0, 13, 0, 21, . . .) , 
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is the sequence a = L^^'^\b) corresponding to the sequence A101890 (see OEIS [9], A000108 , 
Paul Barry 2004). Thus, we know that a = L^~^''^\a) and we find the formulas 



and 




4 Fixed sequences for Interpolated Invert composed 
with Generalized Binomial 

The aim of this section is to find all the sequences fixed by the operator I^^^ o L^^'y\ We 
find a functional equation whose solutions gives us the expression of the ordinary generating 
functions of these fixed sequences and it determines some conditions which the parameters 
x,y and h must satisfy to have solutions. First of all, we recall that for any sequence A G S{R) 
the ordinary generating function L/j(t) for B = L^'^'y\A), satisfies equation (5). The ordinary 
generating function of C = I^^\B) as we defined in equation (1), becomes 

C(t) "^^^'^ 



and we have a fixed sequence if and only if 

A(t) = ^^^^^ ^ ^ (14) 

where A(t) = \A{t) is the ordinary generating function of A and A{t) = Yl o.ni"'^^- Substi- 

n=0 

tuting the expression of L/i(t) in equation (14), with some algebraic manipulations, we find 
the relation 



A 

Ait) 



ht 



h — xA 



ht 
1-yt 



which is equivalent to 



\l-ytj l + xA{t) ^ ' 

Now we study the functional equation (15) remembering that if A G S{R) we have ^4(0) = 
and ^4 (0) = 1. We find explicitly the solutions in the following 

Theorem 7. The operator l'^^^ o L^^'^v) fixes the sequences A G S{R) whose ordinary gener- 
ating functions are of the form 
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if X and h 0,1, and of the form 

A(t) = 
^ ^ 1-ct 

ifx^O,h = l, x + y = 0,ceR. 

Proof. We have to solve the functional equation (15), so differentiating both members of this 
equation with respect to t we obtain 

h f ht \ hA'(t) 



{l-yty \l-ytj {l + xA{t)Y 

multiplying both members by ht"^ and multiplying and dividing the second member by 
{A{t)Y , the previous relation becomes 

ht y ^' f ht \ _t^A'{t)f hA{t) 



1-ytJ \l-ytj {A{t)Y \l + xA{t) 



Now if we use equation (15) to rewrite the second member and we pose w = we find 

..(A(w)? . 

or better 

w^A'iw) t^A\t) 



w'AH='^iAiw)f , 



{Aiw)f iAit)Y ■ 
Finally we come to relation (16) which implies that 

t^A'it) 
iA{t)Y 

where Ci is a real constant. Moreover Ci = 1 because the condition ^'(0) = 1 gives 
Then we can easily solve the differential equation 



(16) 



through separation of variables, obtaining 



Mt) = ^ . (17) 



Substituting the solutions (17) into equation (15) we find 

ht ht 



l-{y + ch)t 1 - (c - x)t 

and this relation is true if and only ii y + ch = c — x oi equivalently c{h — 1) = —{x + y) . 
Clearly if /i 7^ 1 we have c = — frf' if /i = 1 this equality is verified for all c G M if and 
only if X + ?/ = 0. □ 
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Remark 8. We observe that the sequences a which have generating functions satisfying the 
hypotheses of the previous Theorem have the following terms 



or a„, = c" 

respectively when x ^ and /i 7^ 0, 1 or when x ^ , h = 1 , x + y = , c E 



5 Fixed sequences for Revert composed with General- 
ized Binomial and for Revert composed with Inter- 
polated Invert 

In order to determine the fixed sequences for the operator rj o we need to find the 

functional equation related to their ordinary generating function. We recall that the rj 
operator acts as we have described in equations (3), so for a fixed sequence A G S{R) we 
have 

u = tLh{t) , . 

t = uA{u) ^ ' 

where L/i(t) is the ordinary generating function of the sequence L^'^'^^A) defined in equation 
(5) and A(w) = -A{u) is the ordinary generating function of A. Combining the two relations 
(18) we finally come to the desired functional equation 

(19) 



h \l-yt 

The study of this equation and the relative conditions on parameters y and h for the reso- 
lution are summarized in the following theorem 

Theorem 9. The operator rj o L^^^y^ has fixed sequences if and only if h = 1 and their 
ordinary generating functions are 

A{t) ^ 



1 + ft • 



Proof. If we pose in equation (19) w = we obtain the equivalent form 

A (]-Aiw)] = . (20) 

\h ^ ') h + yw ^ ' 

Multiplying both members with |- and applying A to the two sides of this equation we have 
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Now if we use relation (20) with j-A{w) instead of w we find 



and with a httle bit of calculation equation (21) becomes 



A ( ^ ^ = (23) 
\h{h + yw)J h + yA{w) 



Differentiating this relation with respect to w gives 

, , w \ hA (w) 



{h + ywy^ \h{h + yw)J {h + yA{w)f 

and evaluating this equality for w = 0, remembering that A(0) = and ^'(0) = 1, we have 

1 _ 1 

so there are solutions if and only if /i = 1. With this information used in equation (23) we 
retrieve a functional equation 

\l + yw) l + yA{w) 
similar to equation (15) when h = 1, which gives the solutions 

A{w) = . 

1 — cw 

Between these solutions we have to determine the real constant c in order to satisfy equation 
(20) when h = 1. K simple substitution yelds 

A{A{w)) 



1 — 2cw 1 + yw 

which is true if and only if c = — | giving A[w) = ^^y^ - □ 

An analogous result is true when we consider the problem of finding fixed sequences 
A G S{R) for T] o /(^). The functional equation which arises in this case is 

a( ^<'> (24) 



1 - xA{t) 

and comes out combining the following relations 

fA(t) 



U 



l-xtA{t) 

t = uA{u) 

where we made use of equations (3), (1) and A(t) = jA{t) is the ordinary generating function 
of A. All the results concerned with the solutions of equation (24) and fixed sequences are 
resumed in the next theorem. 
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Theorem 10. The fixed sequences for the operator rjo I^^^ have ordinary generating functions 

1 



1 + ft 

Proof. Setting w = ^^^^''^^^ in equation (24), clearly A{t) = j^^, and we obtain 

A{w) = t . 

Applying A to both members of this equation we find 

A{A{w)) = A{t) = 

which is of the same kind of equation (20), that we have already solved , when /i = 1, during 
the proof of the previous theorem. □ 

Remark 11. We observe that the sequences a which have generating functions satisfying the 
hypotheses of the previous Theorem 9 and Theorem 10 respectively have the following terms 

■irif) and a^={-ir('^ 
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